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Abstract 

We elaborate on the color-kinematics duality for off-shell diagrams in gauge theories coupled to matter, by investigating 
the scattering process gg -A ss,qq,gg, and show that the Jacobi relations for the kinematic numerators of off-shell 
diagrams, built with Feynman rules in axial gauge, reduce to a color-kinematics violating term due to the contributions 
of sub-graphs only. Such anomaly vanishes when the four particles connected by the Jacobi relation are on their mass 
shell with vanishing squared momenta, being either external or cut particles, where the validity of the color-kinematics 
duality is recovered. We discuss the role of the off-shell decomposition in the direct construction of higher-multiplicity 
numerators satisfying color-kinematics identity in four as well as in d dimensions, for the latter employing the Four 
Dimensional Formalism variant of the Four Dimensional Helicity scheme. We provide explicit examples for the QCD 
process gg -A qqg. 
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1. Introduction 

Tree-level amplitudes in gauge theories are found to ad¬ 
mit a color-kinematics (C/K) dual representation in terms 
of diagrams involving only cubic vertices, where the kine¬ 
matic part of the numerators obey Jacobi identities and 
anti-symmetry relations similar to the ones holding for the 
corresponding structure constants of the Lie algebra [1, 2], 
as depicted in Fig.l. 

4—^—3 ^3 4^^^^—3 

Figure 1: The Jacobi combination: it can be applied both to Lie 
group structure constants and to the kinematic part of numerators. 

While first studies dealt with the C/K duality within 
scattering amplitudes involving massless partons, more re¬ 
cent investigations pointed to the possibility that such a 
symmetry can be present also when massive particles are 
involved [3-5]. 

The C/K duality, on the one hand, implies the exis¬ 
tence of relations between color ordered tree amplitudes 
in gauge theories and, on the other hand, yields a gauge- 
gravity dual representation of gravity amplitudes, accord¬ 
ing to which they can be expressed as Yang-Mills ampli¬ 
tudes where the gauge-group structure constants are re- 
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placed by a second copy of the color-kinematics dual nu¬ 
merator, [6-35]. When considering multi-loop amplitudes, 
the C/K duality allows to establish relations between the 
(numerators of the) integrands of planar and non-planar 
diagrams. Therefore, it turns into an efficient algorithm for 
generating either high-multiplicity tree-level amplitudes or 
multi-loop integrands, with a better control of the factorial 
growth of the diagram complexity. 

The C/K dual representation was proven at tree-level 
by employing both string theory methods [9, 10, 14-17] 
and on-shell methods [20-22], and it was conjectured to 
hold at higher orders [2]. 

Finding a systematic algorithm to determine C/K-dual 
numerators is not an easy task, because of the wide range 
of transformations, referred to as generalized gauge invari¬ 
ance, underpinning their representation [7, 13, 25, 26, 36]. 
In fact, the search for a C/K-dual representations can be 
formulated algebraically, at least for tree-level numerators, 
in terms of an inverse linear problem, as recently discussed 
in [37, 38]. Nevertheless, it was possible to build an effec¬ 
tive Lagrangian with the property of generating C/K dual 
numerators for tree-amplitudes [39], and non-trivial ex¬ 
amples of dual representation of higher-order numerators 
were found, up to two loops in non-supersymmetric the¬ 
ories, [8, 40, 41], and up to four loops in supersymmetric 
ones, [2, 42-49]. 

In this letter, we study the role of color-kinematics du¬ 
ality within off-shell currents, which enter the construction 
of both higher-multiplicity tree-level and multi-loop ampli- 
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Figure 2: Embedding of the Jacobi combination into either higher- 
point or multi-loop diagrams. 
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Figure 3: Jacobi combination for gg - 
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tudes. We investigate, in a purely diagrammatic approach, 
the origin of possible deviations from the C/K-dual behav¬ 
ior, providing concrete evidence of their relation to contact 
interactions, which was already pointed out to in [1, 50]. 

First, we consider the tree-level diagrams for gg —^ 
W, for massless final state particles, with X = ss,qq,gg^ 
in four dimensions. We work in axial gauge, describing 
scalars in the adjoint representation, while fermions in the 
fundamental one. We deal with the Jacobi relation of the 
kinematic numerators keeping the partons off-shell. Due to 
the off-shellness of the external particles, the C/K-duality 
is broken, and an anomalous term emerges. This anomaly 
vanishes in the on-shell massless limit, as it should, recov¬ 
ering the exact C/K-duality. 

Later, we show that when the Jacobi combination of 
numerators is immersed into a richer topology, associated 
either to higher-point tree graphs or to loop integrands, as 
depicted in Fig. 2, the anomaly corresponds to the contri¬ 
bution of subdiagrams, obtained by pinching the external 
lines of the Jacobi combination. In other words, the Ja¬ 
cobi relation for the numerators in axial gauge, which, in 
the case of on-shell tree-amplitudes, is identically zero and 
resolves the C/K-duality, in the off-shell case, can be ex¬ 
pressed in terms of contact interactions, which we explic¬ 
itly identify for the first time for the processes at hand and 
represent one of the main result of this communication. 
This decomposition, developed in the canonical formal¬ 
ism of Feynman diagrams, shows that the C/K-duality for 
high-multiplicity diagrams naturally holds when the four 
particles entering the Jacobi combination are cut, since, 
in this case, the contribution of the subdiagrams trivially 
vanishes. 

We discuss how our result, which provides a precise 
identification of the anomalies which should be absorbed 
into the redefinition of the trivalent numerators, can be 
used, together with generalized gauge transformations [1, 
50], in order to re-shuffle contact terms between diagrams 
and build on-shell C/K-dual representations for higher- 
point tree-level amplitudes. As successful check of this 
recursive construction, we present the explicit calculation 
for the tree-level contribution to gg —> qqg. 

Finally, we extend the C/K-duality to dimensionally 
regulated tree-level amplitudes. They are the basic build¬ 
ing blocks for the determination of scattering amplitudes 
beyond tree-level within generalized unitarity based meth¬ 
ods, which require trees depending on the regulating pa¬ 
rameter. 

We adopt a novel variant of the four-dimensional helic- 
ity (FDH) scheme [51-53], the so-called four dimensional 


formulation (FDF), recently proposed by some of the au¬ 
thors in [54]. FDF has the advantages of employing a 
purely four-dimensional representation of the additional 
degrees of freedom which naturally enters when the space- 
time dimensions are continued beyond four. We derive the 
C/K relation for the basic four-point currents, and deter¬ 
mine the dual numerators for the process gg —>■ qqg where 
the initial gluons live in d dimensions. 

Algebraic manipulations and numerical evaluations have 
been carried out by using the MATHEMATICA packages 
FeynCalc [55] and S@M [56]. 


2. Color-kinematics duality for scalars 


The process gg —>■ ss gets contributions from four tree- 
level diagrams, three of which contain cubic interactions, 
due to either ggg or gss couplings, while one is given by 
the quartic vertex ggss. Their color factors, for which we 
adopt the normalization = Tr([r“, r^]r=) = 
satisfy the Jacobi identity 
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A similar relation can be established for the kinematic part 
of the numerators of suitably defined graphs involving only 
cubic vertices. 

In fact, after performing the color decomposition and 
some algebraic manipulations, the contribution of the four- 
point vertex can be distributed to the numerators of the 
diagrams with cubic vertices only, hence yielding the iden¬ 
tification of three color-kinematics dual diagrams [1] . The 
corresponding numerators, say ni, n 2 and ns, can be com¬ 
bined in Jacobi-like fashion. 


Ns = -ni + n2+ ns, 


( 2 ) 


as shown in Fig. 3. 

In axial gauge - that we will consider throughout our 
calculations - the numerator of the gluon propagator takes 
the form 

n^''{p,q) = n^:^ + nZiP,q\ (3) 

where Ilp^y corresponds to the numerator of the propaga¬ 
tor in Feynman gauge and g) labels the term de¬ 

pending on an arbitrary light-like reference momentum g^. 
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The explicit form of (2) is given by the contraction of an 
off-shell current with gluon polarizations as^, 

Ns = (j^^Fey + (Pl) ^^4 (Pi) 

= Ns-Fey + Ns-Ax, (6) 

where Feynman gauge-like terms 

of the three numerators, 

-'iJs-FeyiPl^P2,P3,P4) =pT (pf + ‘^P't) 

-(2pf+pf)pf (7) 

and 


- *'^s'^Ax‘‘(P1>P2jP3,P4) = 

- (Pi'Pi^ - Pf P^ - (p? - P 4 ) 9 ■ (P 2 - Pa) 

+ (P 2 - P 3 ) [(P4 + +q - (P4 - Pi) 9^^^^ 

-(pi+2p4)'^^g'^-]| (8) 

is the contribution, depending on the reference momen¬ 
tum, which only originates from 712 . 

These expressions, obtained using only momentum con¬ 
servation X)i=i Pi = Oj show that the Jacobi identity holds 
also on the kinematic side, i.e. Ng = 0, once we im¬ 
pose the on-shell conditions of the four external particles, 
pf = 0, as well as the transversality condition for gluons, 
Pi-e(Pi)=0, 7 = 1,4. 

We want to remark that W-Fey and Ns-Ax vanish sepa¬ 
rately, so that the C/K duality is satisfied at tree-level 
also in ordinary Feynman gauge. 

A similar calculation was performed in [57], where tree- 
level numerators for gg ^ X were studied as well. Our 
result differs in the choice of axial gauge, which, as we are 
going to show, plays an important role in the identifica¬ 
tion of the C/K-duality violating terms in the numerator 
higher-multiplicity graphs. 

The expressions of the currents in Eqs. (7,8) are valid 
for off-shell kinematics. Therefore, they can be exploited 
for providing a better understanding of C/K-duality within 
more complex numerators obtained by embedding the Jacobi- 
like combination of tree-level numerators into a generic di¬ 
agram, as depicted in Fig. 2, where the double circle shall 
represent an arbitrary number of loops and external legs. 

In the most general case, the legs Pi,P 2 jP 3 and p 4 become 
internal lines and polarizations associated to the particles 
are replaced by the numerator of their propagators, which, 
for the scalar case, simply corresponds to a factor i. Ac¬ 
cordingly, Eq. (2) generalizes to the following contraction. 

Ns = {NsUo,X‘^^^\ (9) 


9 ■ (Pi + P 4 ) 


^From now on, the coupling constant as well as any constant 
prefactors associated to the normalization of the generators of the 
gauge group are understood. 



Figure 4: Off-shell color-kinematics duality for gluons and scalars. 
The Jacobi combination of tree-level numerators (l.h.s) is expressed 
in terms of subdiagrams only (r.h.s.). 


between the tensor {Ns)a-iai, defined as, 

{Ns)aia4 = ~{j!^Fey Js-]tx)^ii-iai{Pl, qi) ^fi4a4{P4, 94) ■ 

( 10 ) 

and the arbitrary tensor standing for the residual 

kinematic dependence of the diagrams, associated to either 
higher-point tree-level or to multi-loop topologies. 

Using momentum conservation, we find that the r.h.s. 
of (10) can be cast in the following suggestive form, 

(A^s)qiQ4 = Pl{^s)aia4 + P4{Ag)ctia4 

J" P2('®s)ai«4 J" Pz^^ s)aia4 

+ PlP2(C'l^)aia4 +P?P3(C'l^)«ia4 
+ PlpliCl^)0,40^4 +P4P3(C's'‘)aia4 , (11) 

where A \, Bl and Cl^ are tensors depending both on the 
momenta pi of gluons and scalars, eventually depending 
on the loop variables, and on the reference momenta qi of 
each gluon propagators. 

Remarkably, Eq.(ll) shows the full decomposition of a 
generic numerator built from the Jacobi relation in terms 
of squared momenta of the particles entering the Jacobi 
combination defined in Fig. 3. In particular, this result 
implies that the C/K duality is certainly satisfied when 
imposing the on-shell cut-conditions p| = 0. 

A diagrammatic representation of the consequences of 
the decomposition (11) in (9) is given in Fig. 4, where the 
eight terms appearing in r.h.s. of (11) generate subdia¬ 
grams, obtained by pinching one or two denominators. In 
these subdiagrams A*, Bl and C*! play the role of effective 
vertices contracted with the tensor 

The existence of contact terms responsible for the vi¬ 
olation of the C/K-duality was conceptually pointed out 
already in [1] . Here, we identified, for the first time to our 
knowledge, on a purely diagrammatic basis, the sources of 
such anomalous term, exposed in the (single and double) 
momentum-square dependance of formula (11). 

The choice of axial gauge turned out to be crucial within 
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our derivation, since the p^-terms appear, beside from 
the trivial contraction also from the contraction of 

(p, q) with the corresponding gluon momentum (Ward 
identity), 


p^n^'"(p,g) = ipf,(^-g^'' 


pt^q’'+ qi^p’'\ q^ 

- ] =^P -■ 

q-p ) q-p 


( 12 ) 

For the sake of simplicity, we do not provide the ex¬ 
plicit expressions for Bl, and (7*1. 

By inspection of (7) and (8), we observe that Js-Fey only 
gives contribution to Ag and Af, while Js-Ax produces 
terms proportional to the momenta of all the four particles 
as well as to all the possible pairs of gluon-scalar denomi¬ 
nators, i.e. contributes to all the eight effective vertices. 
In addition, because of the explicit symmetries of Js-Fey 
and Js-Ax under 1 o 4 and 2 o 3, the two effective vertices 
associated to the pinch of one scalar propagator, namely 
Ag and Af, are related to each other by particle relabelling. 
The same happens for and B ^, which correspond to the 
pinch of one gluon propagator. For the same reason, there 
is only one independent Cs function, corresponding to the 
pinch of two denominators, say PiP'j, which are originated 
from terms proportional to in Eq.(8). 


In the following Sections, we show that off-shell color- 
kinematics identities can be established, along the same 
lines, for the coupling of gluons to quarks as well as for 
pure gauge interactions. 


3. Color-kinematics duality for quarks 

The tree-level scattering gg qq has a simpler di¬ 
agrammatic structure than the previously discussed case, 
because of the absence of any four-particle coupling. There 
are three Feynman graphs contributing to it, and they con¬ 
tain only cubic interactions due to ggg and gqq couplings. 
The corresponding color factors obey the Jacobi identity, 

rpa^rpdl I rpairpa4 . ra4aibrT->b_ _ n M 



Figure 6: Off-shell color-kinematics duality for gluons and quarks. 


Using Dirac algebra and momentum conservation, Jq.pgy 
and Jq.A^x tie organized into compact forms as. 


-*'^q-Fey(Pli7’2,P3,P4) = 

= + if 3 + 1 * 2 ) 

+ pf7^*-^^7^*^ (15) 


and 


- * Jq-A^ (pi) P2, P3, P4 ) = 

- q^^ (pi + 2pir^ + g'** (P4 + 2 ^ 1 )^*"] 

- ^ [Pi^pf - pTpT + ipl - Pi)] I • ( 16 ) 

We observe that (14) vanishes when the four external 
particles are on-shell, due to transversality conditions and 
Dirac equation, u{p'i)pf 3 = Pf 2 v{p 2 ) = 0. 

The C/K duality is satisfied at tree-level also in Feynman 
gauge, since on-shellness enforces fVq-Fey and N^-Ak to van¬ 
ish independently. 


and, in this case, it is straightforward to build the com¬ 
bination of color-kinematics numerators for the tree-level 
graph (2), as shown in Fig. 5. 



Figure 5: Jacobi combination for gg -A qq. 


Following the same derivation as for gluons and scalars, 
the Jacobi relation for gluons and fermions can be built 
from the contraction of fermion currents and polarizations, 

Nq = M(p3)(Jq-'Fey + "^q-Ax )^(P2)e,X4 (P4)eAU (Pl) 

= fVq-Fey + fVq-Ax- (14) 


In order to study the Jacobi combination within higher- 
point numerators or multi-loop integrands, we repeat the 
procedure adopted in Section 2. Accordingly, we promote 
the external states of gluons and quarks to propagating 
particles, and define the off-shell tensor, 

(iVq)a4«4 = ¥3«-Ve; + <-Ax )¥2 

X n^4c<4(P4,g4)n^iai(pi,gi), (17) 

by replacing the polarization vectors with the numerators 
of the gluon propagators, and the spinors u{p 3 ) and v{p 2 ) 
with the numerators of fermionic propagators. As be¬ 
fore, the full numerator is obtained contracting (17) with 
a generic tensor 


4 




Manipulating the r.h.s. of (17), we obtain an expres¬ 
sion analogous to the decomposition (11), where the de¬ 
nominators of the four particles are manifestly factored 
out. 


{^q)aia4 —Pl{-^q)aia4 + P oiioti 

+ P2i^q)aia4 d" Pi{^q)aia4 

+ p\pl{CfUo.4+Pipl{Cf)<.4<.4 
+ Plpl{Cf)aia4 + Plpl{Cf)aia4 ■ (18) 

In the above expression, and S* receive contribution 
both from Jq-Fey and Jq_Ax while ’s are determined only 
dy J q-Ax ■ 

This can be understood by inspection of (15) and (16), 
observing that denominators may appear because of (12), 
as well as because of the identity = 

Also in this case, we only have three independent func¬ 
tions: two for the effective vertices corresponding to the 
pinch of one quark- or one gluon- propagator, namely 
and i?*, and a single vertex (7*1 for the pinch of a quark- 
gluon pair. We remark that these functions contain non¬ 
trivial Dirac structure. 

The interpretation of (18) is similar to the one of (11) and 
it is illustrated in Fig. 6. 


-pY 


{P2 + 2paY" 

_ ^Ml/^4 

iP2 

-f 

^paY^ 



+ 9^^^^ {P2 + 2P3Y^ 

1 



+ pY 

|'^M2M4 

(P3 + 2p4)^^ 


{P3 

-k 

2paY^ 



+ g^^^YP3 + ^P2Y" 

1 



-pY 


{pA + 2p3Y" 


{pi 

-k 

2P3Y^ 


Ag^^^^{pA + ^P2T\ (20) 


- i 




{PI,P2,P3,Pa) = 


1 


VAX 

(pfpf - p'i'^p'i^ - {pI - pI) [q • (pa - Ps) 

- {P3 + ‘2‘PaT'' + (pa + 

b - pTpT - {pI - Pa) [q ■ {pi - P 2 ) 9^^^ 


+ {pi + 2p2Y^ q^"" - ip2 + 2 ^ 1 )^"* g'**]I 

- [(1234) ^ (4123)] - [(1234) (4231)]. (21) 


With the by-now usual arguments, we observe that the 
tree-level C/K duality, A^g-Fey = A^g_Ax = 0, holds when 
the external particles are on-shell, separately for the Feynman- 
and axial-gauge contributions. 


4. Color-kinematics duality for gluons 

Finally, we consider the C/K duality for the pure gauge 
interaction process gg —>■ gg. 

As for gg ^ ss, there are four diagrams to be considered, 
three involving the tri-gluon interaction, and one contain¬ 
ing the four-gluon vertex. The color factors obey the Ja¬ 
cobi identity (1). After distributing the contribution of 
the four-gluon vertex into the three structures according 
to the color decomposition, we can define three graphs 
with cubic vertices only whose numerators enter a Jacobi 
combination, as shown in Fig. 7. 




Figure 7: Jacobi combination for gluons. 


Figure 8: Off-shell color-kinematics duality for gluons. 


With this prescription, the kinematic Jacobi identity 
takes the form 

Ng, = {J^Fey^ + '^g-Ay*)eA*l(7’l)eM2(P2)V3(P3)eM4(7'4) 

= A^g-Fey + A^g-Ax, (19) 

where 

- * JY/eT""^Y pi,P2,P3,Pa) = 

= pY {pi + 2 paY^ - {Pi + ‘2 ‘PaY'' 

+ 9^'^^'' {Pi + “^PsYl 


As in the previous Sections, we can build a generic off- 
shell tensor, to be embedded in a more complex topology, 
either with more loops or more legs, by replacing polar¬ 
ization vectors with the numerator of the corresponding 
propagators, 

11 ^ 202 (^ 2 ; 92 ) n^3a3(p3i gs) n^.ja.j(p4, g4) , (22) 

and by contracting this expression with an appropriate 
tensor 




Because of the Ward identity (12), (-/Vg)^,^ turns out 
to be decomposed as 

4 4 

(^g)ai...a4= 

i—1 hj—^ 

(23) 

We observe that, differently from the scalar and fermionic 
cases, Jg-Ax produces all the possible combinations of two 
different denominators. This is a consequence of the per¬ 
mutation symmetry of the gauge-dependent part of the 
numerators, which is an exclusive feature of the Jacobi 
identity for pure gauge interactions. The same symmetry 
reduces from three to two the number of independent ef¬ 
fective vertices, corresponding to the pinch of one or two 
gluon propagators. 

The diagrammatic effect of (23) contracted with 
is depicted in Fig. 8, which shows that, as it happened for 
the gg ^ ss, and gg —>■ qq, the Jacobi combination of the 
kinematic numerators for gg —>■ gg with off-shell particles 
always reduces to subdiagrams. 

Let us, finally, remark that the form factors A®-, B®- 
and C^-^-type appearing in the decompositions (11), (18), 
and (23) still depend on the momenta pi and pj. Therefore, 
within multiloop integrands, they can generate tensor in¬ 
tegrals which can be subject to further integral reductions. 

5. Construction of dual numerators for higher-point 
amplitudes 


whose solution allows us to express M color factors in 
terms oi N — M independent ones {Cc,(i), ..., so 

that (24) can be organized as^ 

N-M 

-4m“(Pl>4'2,---,Pm) = X] (26) 

where 

N 

a.w, e {0,±1}. (27) 

i=i ^ 

The set of identities (25) is not, in general, trivially satis¬ 
fied by the corresponding kinematic numerators, whose Ja¬ 
cobi combinations produce non-vanishing anomalous terms, 

rij J- Tlj J- Tlf^ — 4^[i,j^k]- (^^) 

We observe that the two sets of equations (27) and (28) 
can be conveniently organized into the matrix equation, 

An = K + 4>, (29) 

with 

n = (ni,n 2 , ...,nAr)^, 

K = ({if0,0,...,0)^, 

c/) = (0,0,...,0,{^[„-fe]})^ (30) 

and 

A,, e {0,±l,±Bj-i}. (31) 


In this Section we illustrate how the previous results 
can be used, together with generalized gauge invariance, [1, 
2, 7, 50], in order to explicitly determine dual representa¬ 
tion of higher-point amplitudes starting from Feynman di¬ 
agrams. In addition, we show that our construction allows 
a purely diagrammatic derivation of monodromy relations 
for amplitudes, [10]. 

Any tree-level m-point amplitude can be decomposed in 
terms of cubic diagrams 


N 

1 

1—1 


(24) 


where Ci is the color factor associated to the i*^-graph, Di 
collects the denominators of all internal propagators and 
Ui is the kinematic numerator which, besides the appropri¬ 
ate Feynman rule-term, might contain contributions from 
contact interactions, that are assigned with the prescrip¬ 
tion described in Sections 2 and 4. 


As shown by the decompositions (11), (18) and (23), the 
anomalies are proportional to the off-shell momenta of the 
particle entering the Jacobi combination itself. Therefore, 
the rise of these anomalies seems to be related to the al¬ 
location of contact terms between cubic diagrams, which 
naturally provides numerators satisfying C/K-duality in 
the four-point case only. 

As a consequence, in order to obtain a dual representa¬ 
tion of the amplitude, we need to re-shuffle contact terms, 
leaving (24) unchanged. This can be achieved through a 
generalized gauge transformation, which consists in a set 
of shifts of the kinematic numerators, 

n* -)■ n' = Ui - Ai, (32) 


satisfying 


N-M N 

(PljP2j '‘',Pm) = ^ ^ ^a{i) ^ ^ 

i=l j=l 


0:cr(z)j Aj 

D, 


0, (33) 


The N color factors appearing in (24), satisfy a set of 
M, M < N , Jacobi identities 

—Ci Cj -|- Cfc = 0, (25) 


^ The amplitude decomposition in terms of a color basis obtained 
as a direct solution of the Jacobi/Lie algebra relations is equivalent, 
yet different from the decomposition proposed in [58]. In particu¬ 
lar, the kinematic coefficient of each color factor appearing 

in (26) is gauge invariant, and corresponds, in general, to a linear 
combination of color ordered partial amplitudes. 
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in such a way that the amplitude can still be written as 

N , N-M 

2=1 2=1 

(34) 

By imposing the vanishing of the coefficient of each CCT(i) in 
(33), the gauge invariance requirement is translated into a 
set equations for the shifts, which leaves M of them unde¬ 
termined. This means that, in principle, we have enough 
freedom to ask the shifts to be solution of M additional 
equations, 

— Ai -I- Aj -I- Afc = (35) 

which, inserted in (28), make the new set of numerators 
n' manifestly dual. 

Thus, the simultaneous imposition of (33) and (35) leads 
the determination of numerators satisfying the C/K-duality 
back to the solution of the N x N linear system 

AA = 0, (36) 

with 

A= (Ai,A2 ,...,A^)^, (37) 

whereas the vector (p and the matrix A are the ones de¬ 
fined by (30) and (31). 

By solving (36), we can determine the shifts to be per¬ 
formed on the numerators obtained from Feynman dia¬ 
grams ensuring C/K-duality as a function of anomalies 
4>[i,j,k] and denominators Di. 

We note that the existence of a dual representation of 
the amplitude is bound to the consistency of the non- 
homogenous system (36), i.e. to the condition 

rank(A|^) = rank(A), (38) 

where A10 is the augmented matrix associated to A. 

In particular, if the system had maximum rank N, the ex¬ 
pression of the numerators would be completely fixed by 
C/K-duality. 

However, as we will show in an explicit example, the rank 
of the system turns out to be smaller than N, so that its 
solution will depend on a set of arbitrary shifts, which are 
left completely undetermined by the imposition of C/K- 
duality. The existence of a residual freedom in the choice 
of the dual representation was first observed in [1] and 
more recently, in [37], where the reduction of the tree- 
level C/K-duality to an underconstrained linear problem 
is addressed in terms of a pseudo-inverse operation, it has 
been interpreted as the hint of a possible analogous con¬ 
struction at loop-level. 

We observe that, if 

rank(A) < A, (39) 


the condition (38) can be satisfied only if a number N — 
rank(A) of relations can be established between the anoma¬ 
lous terms ■ 

In the following Section, we will show that these con¬ 
straints, which were obtained in [10] for the five-gluon am¬ 
plitude using string-derived monodromy relations, are fully 
implied by the diagrammatic expansion of the amplitude 
and that they can be obtained from the simple knowledge 
of the matrix A. In particular, they are found by deter¬ 
mining a complete set of vanishing linear combinations of 
rows of A. Working on a specific example we will argue 
that our off-shell decomposition (11),(18) and (23) make 
these relations manifest. 

Moreover, we observe that, because of (29) and (36), ap¬ 
plying the matrix A to the new set of numerators n(, we 
obtain 

An' = K, (40) 

with K given by (30). 

Again,if rank(A) < N, the consistency condition 

rank(AjK) = rank(A), (41) 

implies the existence of A — rank(A) constraints between 
the kinematic factors which are in one-by-one cor¬ 

respondence with the relations between color ordered am¬ 
plitudes, that have been previously conjectured as an im¬ 
plication of of C/K-duality [1] , and then derived from the 
low energy limit of string theory [10] and on-shell recur¬ 
sion [20]. 

Therefore, this construction shows that all these non-trivial 
relations can be derived from the expansion of the ampli¬ 
tudes in terms of Feynman graphs through purely algebraic 
manipulation on the matrix A. 

Finally, we want to remark that, whereas in the usual 
top-down approach, the numerators appearing in the r.h.s. 
of (24) are interpreted as abstract re-organization of Feyn¬ 
man rules-numerators on which, by assumptions, the C/K- 
duality is imposed, our approach provides a systematic 
way to identify the link between dual numerators and 
Feynman diagrams. Starting from the set of explicitly 
C/K-violating but well-defined Feynman rule-numerators 
we determine, by mean of generalized gauge transforma¬ 
tions, the actual redistribution of contact terms which has 
to be performed in order to establish the duality. 

In this framework, the off-shell decomposition of the four- 
point identities derived in the previous sections plays a 
key role in the identification and in the algorithmic con¬ 
struction of the anomalous terms, i.e. the non-vanishing 
element of the vector 0. In fact, the l.h.s. of each kine¬ 
matic equation (28) can be obtained from the contrac¬ 
tion of {Nk)ai...aj {k = s,q,g, depending on the process 
under consideration) evaluated on a suitable permutation 
of the labels of the external legs, with lower-point func¬ 
tions. Therefore, all the anomalies Plij.k] can be deter¬ 
mined, without going trough the explicit calculation of all 
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N numerators, just by identifying the M tree-level subdi¬ 
agrams that can be factored in each of three numerators 
appearing in (28). 

Summarizing the diagrammatic approach to the construc¬ 
tion of C/K-dual numerators for higher-point amplitudes: 

- given the decomposition of an amplitude in terms 
of Feynman diagrams, it is organized into N cubic 
graphs, whose numerators satisfy the system of equa¬ 
tions 


An = K-f0. (42) 

- A generalized gauge transformation 

u-i —^ -f Ai, (43) 

such that 

AA = cj), (44) 

An' = K, (45) 

is performed on the amplitude in order to obtain 
a new set of numerators satisfying the C/K-duality. 
The solution of (44) determines the shifts linking the 
starting set of numerators to the dual representation. 

- The existence of solutions for the the systems (44) 
and (45) is related to the constraint 

rank(A|0) = rank(A|K) = rank(A). (46) 

This consistency condition is able to detect all N — 
rank(A) non-trivial constraints both between the C/K- 
violating terms 4’[i,j,k] ^^nd the kinematic factors , 
the latter corresponding to the well-known relations 
between color ordered amplitudes which were first 
observed, for gluon amplitudes, in [IJ.Note the N — 
rank(A) also determines the number of completely 
free parameters the set of C/K-dual numerators will 
depend on. 

In the following Section we give an example of this method, 
determining the C/K-dual representation for gg —>■ qqg 
and showing that the knowledge of the matrix A can be 
used to determine the constraints on kinematic factors 
Ka(i) ss well on the anomalies i which rise as direct 
consequence of the off-shell decompositions worked out in 
Sections 2-4. 





Figure 9: Feynman diagrams for gg —>■ qqg 


The process under consideration contains a single exter¬ 
nal quark-antiquark pair and, as a consequence, receives 
contribution from the four-gluon vertex. This allows us 
to show, in a concrete case, how contact interactions can 
be treated. We go step by step through the procedure 
outlined in Section 5, adopting notation and conventions 
similar to [1]. 

Fig. 9 shows the 16 Feynman diagrams for the process 
gg —>■ qqg. The contribution of nig, which, containing the 
four-gluon vertex, depends on three different color struc¬ 
tures, 


cie’^ie = Canape + c^n^-ie + cangpe (47) 

can be split between na, ng and ng, so that the new cubic 
numerators read 


na + Si2na;i6 -)■ na, 
ng + sisngpe —>■ ng, 

ns + S25''^3;16 ^s, (48) 

being Sy = {pi +Pj)‘^. 

Thus, the decomposition of the amplitude in terms of cubic 
graphs reads 


6. Color-kinematics duality for gg —>■ qqg 

Extensions of C/K duality in QCD amplitudes with 
fundamental matter have been discussed in [3, 5], where 
manifest duality has been verified for several processes. 

In order to illustrate the method proposed in the previous 
Section, we provide a further example of C/K-duality in 
QCD, by determining dual numerators for gg —>■ qqg. 


4-'^(5(l),5(2),9(3),g-(4),5(5)) = 


cini 


C2n2 


C3«3 


C4n4 


Ceng 


S12S45 S32S51 S34S12 S45S23 S15S34 

^ ceue ^ CTTi’j ^ Cgng ^ cgug ^ cipnip 
S14S25 S23S14 S25S34 SiaS25 S24Sia 

_l_ ciinii ^ Ci2ni2 ^ ciania ^ Ci4ni ^ cignig 
S15S24 Sl2Sa5 Sa5S24 Si4Sa5 SiaS45 ’ 


( 49 ) 



















where only 713, and ng differ for an additional contact 
term from the expression given by Feynman rules. 

The color factors, 


„ _ £aia 2 brpb_rpa 5 

— J -^3k^k4: ’ 

_ _ n^b_ xba^c xcaic 

<-3 — -^34/ / 

_ ^a^aibJba 2 Crpi 


34’ 


= 

c, = T^iTlJ^ 

„ _ £a 5 aibrpb_rpa 2 

(-11 — J 

C 13 =^ 3 “|T-T“^ 

ci 5 = 


^ _ rpa 2 rpb_ fba^ai 

^2 - ^Sk^k^J 

^4 = -r3|r-r;/, 

C6 = T3\T“i/ 

C8 = -r^“= 7 '’ - 34 > 

„ _ rpairpazrpa 2 

^10 “ ^3fe^fc/ ’ 

rjia^rpb fba\a 2 

Cl2-J3fc^fe4/ 


b rj^ai fba 2 (i 5 

7 

0 , 2 Cl 5 b fb(l\CrpC 


3k kl M 

satisfy a set of 10 Jacobi identities, 


(50) 


— Cl + C3 + C12 = 0, 

— Cl + C4 + C15 = 0, 

— C2 + C4 + C7 = 0 , 

— C2 + C5 + Cii = 0 , 

— Cg + C7 + Ci4 = 0 , 

— Cg + Cg + Cg = 0 , 

— Cg + Cio + Ci5 = 0 , 

— Cii + Cio + C13 = 0 , 

— C12 + Ci3 + Ci4 = 0 , 

(-C5 + C3 + Cg = 0 ). ( 51 ) 


The system is redundant, since any of the above equa¬ 
tions, for instance the last one, can be expressed as a lin¬ 
ear combination of the others 9 . Therefore, it can be freely 
dropped. 


We solve (51), choosing {ci,cg, C3,C4, cg,cg} as indepen¬ 
dent color factor, and we re-express the amplitude as 

6 

4 ’--(g(l), 5 ( 2 ),g( 3 ),g( 4 ),ff( 5 )) =^c,JC, ( 52 ) 

i=l 


Kl 

K2 

Ks 


Ki 

K5 

Ke 


ni_ ni2 ni 3 


S12S45 S12S35 S24S35 

712 n-r 7714 

S23S51 S14S32 S14S35 

ns ng _ 7712 

S12S34 S13S25 S12S35 

7^4 Tlj 71 i 4 

S23S45 S14S32 S14S35 

775 77 g 7710 

S34S51 S13S25 S13S42 

776 ^ 779 ^ 7714 

S14S25 S13S25 S14S35 


7710 ^^15 


S13S42 

-I- 

S13S45 

1 ni 3 


7711 

S24S35 

“T 

S42S5I 


ni3 


nio 

ns 

S24S35 


S13S42 

S25S43 

ni3 


nio 

ni5 

S24S35 

“T 

S13S42 

S13S45 

1 77 g 


7711 


S25S43 


S42S5I 

J 

ni3 

-f 

nio 

. (53) 


S24S35 S13S42 


We remark that this choice is by no means unique and 
other admissible sets of independent color factors lead to 
different but equivalent decompositions. 


Now we introduce the shifts ( 32 ) and, by imposing gen¬ 
eralized gauge invariance on ( 52 ), = 0 , we obtain a 

set of 6 homogenous equations. Furthermore, in order to 
establish C/K-duality for the new numerators 77', we re¬ 
quire the shifts to absorb the anomalous terms, 7 .e. to be 
solution of an additional set of 9 non-homogeneous equa¬ 
tions, that are obtained from ( 51 ) by replacing each color 
factors Ci with the corresponding shift and the r.h.s. 
with the proper anomaly. 

Thus, dual numerators are determined from the solution 
of the linear system ( 36 ), where A is the 15 x 15 matrix 


J 


S12S45 

-1 

-1 

0 

0 

0 

0 

0 

0 

0 


S23S15 

0 

0 

0 

0 

0 

-1 

-1 

0 

0 

0 

0 

0 


S12S34 

0 


S34S15 

0 


S23S45 

0 


S14S25 

0 

0 

0 

0 

0 

0 

-1 

-1 

0 

0 

0 


S14S23 

0 

1 

S14S23 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 


S25S34 

0 


Sl^S25 

S13S25 

0 

1 

S13S25 

0 

0 

0 

0 

0 

0 

0 

1 

-1 

0 

0 


■S13S24 

0 

1 


’S13S24 

0 


424S15 

0 

1 

S24S15 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

-1 

0 


S12S35 

0 


S12S35 

1 

0 

0 

0 

0 

0 

0 

0 

-1 


24jS35 


■S24S35 

0 


S14S35 

0 


and the vector cj) is given by 

4 > = (0, 0, 0, 0, 0, 0, 4 >[ i ^ s ^ i 2 \, <('[1,4,15]) <('[2,4,7] 5 ^^[2,5,11]) ^^[6,7,14]) ^^[6,8,9] j ^^[9,10.15] j <('[11,10,13] j ?^[12,13,14])"^- 


S 13 S 45 

0 


(54) 


( 55 ) 





































































The anomalies 4>[i,j,k] can be directly obtained from our 
off-shell decompositions. We observe that, since we de¬ 
cided to drop, without loss of generality, the last of equa¬ 
tions (51), which involves a Jacobi identities for gluons, we 
can express all anomalies in terms of the fermionic current 


dq = Jq-Feyn + Jq-Ax, whose explicit cxpression is given 
by equations (15) and (16). Likewise, the three anomalies 
involving the numerators na, and ng receive an addi¬ 
tional contribution from the four-gluon interaction, since 
no contact term was considered in the definition of Jq, 


^[1,3,12] 

^^[ 1 , 4 . 15 ] 

<(>[2,4.7] 

?^[2.5.11] 

^^[6,7.14] 

<(>[6,8,9] 

<(>[9,10,15] 

<(>[11,10,13] 

0[12,13,14] 


U3 [jq (pi +P2,P4,P3,P5)aa^] {pi -b P2, <?) Vsai as (“Pi - P2, Pi, P 2 ) sT + «3;16 = Sl2‘b>[1.3,12], 


U 3 


U 3 


Jq (P2,P4 +P5,P3,Pl)„2ai k 

k (^ 3 +^ 2 ) '^q(?’l>P4,P3 +P2,P5)o,a5 


V4E°1^£^'^ — 545‘/2[1,4,15] j 


U4£“^£5® = S23</>[2.4,7], 


M3 [^q {pi +P5,P4,P3,P2)aa2] '^4 (pi +P5,q) Vsaias (“Pi - P5, Pi, Pd) sT k"" + n-5;16 = Sl5<P[2.5.11], 


U 3 


Jq iPd,P4+Pl,P3,P2)c,,a^ {h+k) k 


VieTk'" = Sl4<b>[6.7,14], 


U 3 [Jq {P2+P5,P4,P3,Pl)aai] ^4^1“^ {P2+P5,q) Vsasas {-P2 - Pb,P2,Pb) k'" k"" +’^8;16 = S25<Pl6.8,9], 


U 3 

U3 

U3 


k (^3 +^l) Jq{Pd,P4,P3 ^Pl,P2), 

Jq (Pl,P4 +P2,P3,P5)„ia, (^ 4+^2 

k ( 1^3 +(^ 5 ) kk>i,P4,P3^Pd,Pi)^^^^ 


1;4£2»£5® — S13</>[9,10,15] j 
U4£“^£5® = S24<P[11,10,13] I 


'C4£i^£2^ — S35</>[12,13,14]) 


(56) 


where V^-^^^fi^(pi,p 2 ,P 3 ) stands for the kinematic part of 
three-gluon vertex. 

We remark that in (56) the factorization of Mandelstam 
invariants follows from the decomposition (18) and 
that, for the anomalies with an internal gluon propagator, 
namely ?i>[i, 3 p 2 ]! <(>[ 2 . 5 , 11 ] and <(>[ 6 , 8 , 9 ]; it can be achieved in 
a straightforward way thanks to the choice of axial gauge. 

As we have already anticipated, the system of equations 
is redundant. In fact, by using momentum conservation 
to express all the invariants Sij in terms of 5 independent 
ones, for instance {si 2 , S 23 , S 34 , S 45 , S 51 }, we obtain 

rank(A) = 11. (57) 

Therefore, if a solution exists, there must be constraints 
between the non-zero elements of (f> able to lower the rank 
of the adjoint matrix. In particular, we expect these re¬ 
lations to correspond to four independent vanishing linear 
combinations of rows of the matrix A. 

This observation provides a constructive criterion to find 
out the constraints between anomalous terms. 

First, we build the most general linear combination of rows 
of the matrix A and we fix the coefficients by requiring 

15 

^AA, =0. 

i=l 


According to (57), one can find at most four linear inde¬ 
pendent solutions {/3p^} to (58). 

Secondly, after selecting an arbitrary complete set of solu¬ 
tions j = 1, 2,3,4, we can verify that 

15 

^/3pV. = 0, Vj = I,2,3,4, (59) 

i=l 

that gives the desired constraints between the C/K-violating 
terms. 

For this specific case we find. 


<(>[6,7,14] 

1 <(>[9,10,15] 

^[12,13,14] 

^^[1,4.15] 

= 0 

Si4 

S13 

S 35 

S 45 


<(>[6,8,9] 

1 <(>[9.10.15] 

<(>[12,13,14] 

<(>[1,3,12] 

= 0, 

S25 

Sl3 

S 35 

S 12 


^^[6,7,14] 

1 <(>[9,10,15] 

^[11,10,13] 

^^[2,4,7] 

= 0, 

Si4 

Sl3 

S 24 

S 23 


<(>[6,8,9] 

<(>[11,10,13] 

1 <(>[2,5,11] 

<(>[6,7.14] 

= 0, 

S25 

S24 

Sl5 

Si4 



which, thanks to (56), can be written as 

‘7’[6.7.14] + A[9.10,15] ~ <^>[12,13,14] ” </>[!.4.15] = 0, (61) 

‘7’[6.8,9] + <^[9,10,15] ~ <^[12,13,14] “ <P[1.3,12] = 0; (62) 


(58) 
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V[Q,7M] + ‘/’[9,10.15] ~ <^[11,10,13] ~ ‘^[2.4.7] — 0; (63) 

'/7[6,8,9] + '/7[ii40,13] + V7[2.5,11] ~ ‘^[6,7,14] = 0. (64) 

These relations follow as a direct consequence of the off- 


shell decomposition (18). For instance, if we specify it for 
the anomalies appearing in (61) by using the explicit form 


of Jq, we obtain 


V7[1,4,15] = “3 
V7[6,7,14] = M3 
‘^[12,13,14] = “3 
= “3 


_ gO!loi 2 ^as _ 
_ gOi^a^^ai _ 
_ ^«5^0!i^a2 _|_ gaia2^a5 
_ ^011^02^015 _|_ gOt2a5^ai 


gaiO‘ 2 q ■ - 

q ■ {pi +P 2 ) 

gO‘20‘nq ■ _ 2 {q°“^p^^ - 

q ■ {p2+Pb) 

^ g°‘^°‘^q ■ {pi -P 2 ) - 2(g“^Pi^ 

q ■ (pi +P 2 ) 

^ 9°“^°‘^q ■ {P2-Pb) - ‘^{q°“^pT 

q ■ (p2 +Pb) 


n0‘2^0lj_\ 1 

g P 2 ) „,a 5 

05 Ot 2 \ 1 

g P 5 > „,ai 

-q^^pT) 

-7 

-q^^pT) 

-7 


^4'^Q:i^Q'2 ^ 0:5 5 
^4'^Q:i^Q'2 ^ 0:5 5 

^4‘^ai 1^02 j 

^4‘^ai 1^02“-as ■ 


(65) 


By using Clifford algebra we can verify that 

</’[!, 4 ,15] + <P[12,13,14] — P[&,7M\ — </7[9.10,15] = 

2£r„iea2^a5“3[5“"“''7“® - = 0 . ( 66 ) 

Similar cancellation are encountered in all other cases. 
The constraints (60) make the consistency relation satis¬ 
fied, 

rank(A|^) = rank(A) = 11. (67) 

As a consequence, the system admits a solution which 
leaves four shifts completely undetermined and the am¬ 
plitude has a C/K-dual representation, consistent with 
generalized gauge invariance, whose numerators depend 
of four free parameters. This number agrees with the 
[n — 2)! — (n — 3)! degrees of freedom found in [1] for 
the pure Yang-Mills case. 

In order to find an explicit expressions for the shifts, 
we build a maximum-rank system by selecting a subset of 
11 independent equations and proceed by Gaussian elimi¬ 
nation. We observe that, as for the solution of the Jacobi 
identities for color factors, also in this case there is a re¬ 
markably large freedom in the choice of the independent 
equations to be solved and, furthermore, in the set of four 
arbitrary shifts to appear in the solution. 

In our case, by selecting equations corresponding to rows 
1-2 and 9-15 of (54), we express A^, z = 5,6,..., 15 as lin¬ 
ear combination of the anomalies of the four 

arbitrary shifts {Ai, A 2 , A 3 , A 4 }, 

4 

Aj = TT-ij (sfci)Aj-l- TZi 

[m,n,p] i^^kl ) 0 [m,n,p] i 

J —1 [m,n,p] 

for i = 5, ...15, ( 68 ) 

where Ti-ij and TT-q^np] are dimensionless rational functions 
of the invariants su ■ 


The analytic expression of (68) , which is not provided here 
for sake of simplicity, has been obtained for arbitrary po¬ 
larizations and has been numerically checked for all helicity 
configurations. In particular, the complete independence 
on the actual values of the four independent shifts has 
been verified for the full color-dressed amplitude as well as 
for each ordering appearing in (52). We observe that the 
choice Ai = 0, i = 1,...,4 leads to a dual representation 
where four numerators correspond exactly to the starting 
ones and three anomalous terms are attributed to single 
diagrams 

^7 = 4>[2,A,7]i 
Ai2 = <(>[1,3,12]) 

Ai5 = <(>[1,4,15]. (69) 

In addition, for any choice of the free parameters, the set of 
new numerators n', satisfies the system of equations (40), 
where 

K = (iFi,X2,...,i^6,0,0,...,0)^. (70) 

Therefore, the consistency requirement 

rank(A|K) = rank(A) = 11 (71) 

allow us to use the exactly the same solutions of 

(58), to establish relations between the kinematic factors, 

15 

^/3p)K, = 0, Vj = l,2,...,4, (72) 

i=l 

which, in this case, read 

SibKi — S34K3 — S14KQ = 0 , 

5127^1 — S 23 K 4 — S 23 KQ = 0, 

Sl5Ar2 — S 43 K 4 — S 25 KQ = 0, 
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S23K2 — S34K5 + (S 23 + S35)Kq = 0 . ( 73 ) 

This set of constraints reduces from 6 to 2 the number of 
independent kinematic factors. 

Therefore, as we have already pointed out, (73) can be 
considered as equivalent to the well-known monodromy 
relations which have been shown, for the pure-gluon case, 
to reduce to (n — 3)! the number of independent color or¬ 
dered amplitudes. 

Nevertheless, we want to remark that, in the approach we 
have presented, the origin of (73), as well as the one of 
(60), is shown to be purely diagrammatic. In particular, 
whereas in [10] monodromy relations analogous to (73) are 
derived from the field-limit of string theory and a set of 
relations equivalent to (60) is presented as a parametriza- 
tion their solution, here both are derived as a necessary 
consequence of the redundancy of kinematic matrix A and 
they have be shown to naturally emerge from the off-shell 
decomposition of the Jacobi-combination of kinematic nu¬ 
merators in axial gauge. 

7. Color-kinematics duality in d-dimensions 

In this Section we study the C/K-duality for tree-level 
amplitudes in dimensional regularization. They are the ba¬ 
sic building blocks for the determination of higher-order 
scattering amplitudes within generalized unitarity based 
methods. We employ the four dimensional formulation 
(FDF) scheme, recently introduced in [54]. Within FDF, 
the additional degrees of freedom which naturally enters 
when the space-time dimensions are continued beyond four, 
such as spinors and polarizations, admit a purely four¬ 
dimensional representation. FDF has been successfully 
applied to reproduce one-loop corrections to gg —>■ gg, 
qq —^ gg, gg ^ Hg (in the heavy top limit), as well 
as gg —?► gggg [59]. Accordingly, the states propagating 
around the loop are described as four dimensional mas¬ 
sive particles. The four-dimensional degrees of freedom of 
the gauge bosons are carried by massive vector bosons (de¬ 
noted by g*) of mass g (associated to three polarization 
states) and their {d — 4)-dimensional ones by real scalar 
particles (s*) of mass g, being g an extra-dimensional 
mass-like parameter. A d-dimensional fermion of mass m 
is instead traded for a tardyonic Dirac field (q*) with mass 
m -I- ig-j^ [60] (associated to two spinor states). The d 
dimensional algebraic manipulations are replaced by four¬ 
dimensional ones complemented by a set of multiplica¬ 
tive selection rules. The latter are treated as an alge¬ 
bra describing internal symmetries. In Appendix A, for 
completeness, we provide the Feynman rules of the FDF 
scheme. 

We anticipate that the C/K duality obeyed by the nu¬ 
merators of tree-level amplitudes within the FDF scheme 
are non-trivial relations involving the interplay of massless 
and massive particles. 



Figure 10; Feynman diagrams for g*g* —> gg. 


7.1. Tree-level identities in d-dimensions 

As a starting point, we focus on four-point tree-level 
amplitudes, showing that C/K-duality can be established 
for all interactions involving particles propagating within 
the FDF framework. 

Envisaging one-loop applications of these results, we con¬ 
sider amplitudes with two d-dimensional external particles 
and two four-dimensional ones. Following the by now fa¬ 
miliar procedure of Sections 2, 3 and 4, we build an off-shell 
Jacobi-like combination of kinematics numerators for each 
of the seven processes involving, according to the Feynman 
rules of Appendix A, two external massive particles. We 
show that, in every case, C/K-duality holds after physical 
constraints are imposed. 

a) We start from the scattering of two generalized gluons 
producing two massless ones, g*g* —> gg, with on-shell 
conditions = p\ = 0 and P 2 = P 3 = g^- The am¬ 
plitude receives contribution from the four diagrams 
shown in Fig. 10, where dotted lines indicate generalized 
particles. Exactly as in the massless case discussed in 
Section 4, the four-point vertex contributes to the same 
color structures as the trivalent diagrams. Hence it can 
be decomposed as 

04714 = Cini;4-I-C2n2;4 + C3C3;4, (74) 

so that, by absorbing its kinematic part into the nu¬ 
merators of the hrst three diagrams, the amplitude can 
be expressed in terms of cubic graphs only. 

We observe that 712 is associated to the massless de¬ 
nominator (pi paY while Til and 773 sit, respectively, 
on the massive propagators {(jPi-\-p 2 )‘^ — g^^ and ((piJ- 
Psf - d^)- 

Therefore, the proper definition of the cubic numera¬ 
tors absorbing the contact interaction is 

ni -b ((pi -bP2)^ - M^)ni;4 Til, 

n2 + (Pl +P4)^n2;4 ^2, 

n3 + ((Pi + Pa)^ - d^)n3;A n^. (75) 

With this prescription, the Jacobi-combination of the 
kinematic numerators depicted in Eig.ll can be written 
as 

A^g*g*gg = -'pg.gg£ 7*1 (pi)£/x2 {P2)e^„ (P3)eM4 {Pi)^ (76) 
with 


tMiM2A^3A^4 _ tMiM2M3M4 i t/^1M2M3A^4 i tM1M2M3M4 

•^g’g’gg “"'g-Fey "'g'g'gg-Ax “I" -^g.g.gg.^s- ) 
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Figure 11: Jacobi combination for g*g* —>• gg. 


• 1 ^^JULSLSJUc' 

.0 4^^ * ''3 



2 1 
3 4 


ni n2 ns ua 

Figure 12: Feynman diagrams for gg —>■ s‘s‘. 


•2 

3 


In Eq. (77) Jg-Fey is the same as in (20); Jg«g«gg-Ax is 
given by 

-*Jg.g.gg-Ax= 

(pfT’f - Pi'Pi'" + {p\ - pI) [q ■ {p 2 - Pa) 

+ (2p3+P2)''^ q^^-i2p2+pa)^^ q^ 

+(pf-pf - pf + (P 2 - Pa) [i-bi - Pi) 

+ (2pi +P4)^" - (2P4 +Pi)''^ (78) 


where the off-shell current is 

tMiM4 _/ tMiM4 j_ Tgig4\^AB /'QQ^ 

‘^ggs’s* ~ Wg-Fey ^g-Ax 

We observe that, although ni and na depend on the 
mass of the scalar particle, the /i-dependence cancels in 
their combination. Therefore, C/K-duality, i.e iVggs»s* = 
0, holds in this case as well as in the massless one, ac¬ 
cording to Eqs. (7) and (8). 


whereas the //-dependent term reads 


_ j jMlA‘2At3M4 _ 

■ ■ s-g^ 


1 


s's'&n-g /j,^ 

ggigsgg.gA ((^2 _ ^ 2 ^ ^^2 _ ^ 2 ^ _ ^ 4 ^ 

{{pI-pI){pI-pI)-p^) 


, ggig2g^sgi 


- (P4 - 

pI) 

^M2M4 

(pf 

pf- 

pfpf) 

+ (P?- 

pD 


(pf 

p^- 

pfpf) 

-{pI- 

pI) 


(pf 

p^- 

pfpf) 

+ (P4 — 

pI) 

^M3M4 

(pf 

pf - 

pfpf) 

+ (P3^>f -Pfpf 

)(pfpr 

-pfpf 

+ [pApA 

-pApT) 

1 (pf pf - 

-pfpr) 


(79) 


By inspection of (79) and (78) we see that C/K duality, 
which corresponds to IVg.g.gg = 0, is recovered because 
of transversality, Si-pi = 0, and on-shellness conditions. 


b) Now we consider gg —>• s*s* {pf = pi = 0, = Pi = 

p^) amplitude, whose diagrams are shown in Fig. 12. 

In this case, since the four-point interaction only con¬ 
tributes to two color structures, 

04714 = Cini;4 -I- C3C3;4, (80) 

we can absorb its kinematic part in the two diagrams 
involving a massive scalar propagator, through the sub¬ 
stitutions, 

ni -I- ((pi -I-P2)^ - P^)ni;4 Til, 

«3 + ((pi +P3)^ - P^)ll3;4 113, (81) 

whereas 712 stays the same as defined by Feynman rules. 
In this way, the Jacobi combination of the cubic numer¬ 
ators, depicted in Fig. 13, becomes 

iVggs«s» = (pi)e^^ (p4), (82) 



- 2 . 2 




Figure 13: Jacobi combination for gg s*s*. 


The remaining processes, whose tree-level identities are 
depicted in Fig.l4.(c)-(g) do not involve contact interac¬ 
tions, so that the construction of their Jacobi-combinations 
directly follows from the Feynman rules of Appendix A. 

In order to avoid repetitive discussion, we simply list the 
results giving, for each process, the corresponding on-shell 
conditions and the expression of the Jacobi-combination 
in terms of off-shell currents. 

According to the case, C/K-duality is recovered once transver¬ 
sality of the gluon polarizations and Dirac equation are 
taken into account. We recall that, for a generalized gluon 
g* of momentum p, one has Si ■ p = 0 (i = ±, 0), while for 
generalized quarks q* of momentum p, one has u(p){'^ — 

= 0 and {p -|- ip"f^)v{p) = 0. 

c) 5*5* 55 (P? = P4 = pI=pI = 0): 

iVg.g.qq = u{pa)J^."^.\^v{p2)e^, {Pi)£gAPi), (84) 

with 

_ 77*17*4 t7*i7*4 _i_ t7*i 7*4 78f;'l 

'^g*g*qq "^q-Fey ''q-Ax ’ 1°'-'/ 

where the terms in the r.h.s. are defined by Eq. (15) 
and (16). 

d) qq {pj = pj = pI = pj = 0): 

-^s^s'qq — '^{Ps) ^s* s* {p 2 ^ ■; 

with 


iPl+P4:)-q^ 

5-(P4-Pi) (^ 2 +^ 3 ) - (P4-P?) (87) 
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e) gg 9*g* {pI =pI = Q,pI=pI = p^)- 


-^ggq*q* “ 'fJ‘iP3)Jggq^q’'^(j'2), (88) 

with 

- * -^ggcf-V = -{h- 

- + (^3 +f^ + ipi^^ g^^^^ 

+ P^4^r^-P^^r^ + J^!^:- (89) 
f) g'g -t 9*9 {pI =pI= A^^ pI=pI = O): 

TVg.g-.q = U(P3) Jpg'q‘.q7;(p2)£Ml (Pl)^A‘4 (P4), (90) 

with 


- i Jg."gq.q = - 7 ^ 17 '"" (^2 + 

+ (^3 + ^2 + + P4^7^" - P^l^^ 

+ ^ (^3 + ^2 + iPi^Pi" - pTpT) ■ (91) 

g) s'g q'q {pi = pi = / 4 ^ pI=pI = 0 ): 

iVs*gq*q = u(p3)^s''*gq*q^(l^2)e,X4(P4), (92) 

with 

-*-li"*gq*q = (^2 +*^^^)+^'4'‘7®^r^- 

(93) 

In summary, we have shown that C/K-duality can be 
extended to four-point tree-level amplitudes in FDF, which 
would naturally enter the construction of loop-level ampli¬ 
tudes in the framework of d-dimensional unitarity. This 
constitute the second main result of this communication. 

In the concluding part of this letter, as a further applica¬ 
tion of the method described in Section 5, we will derive 
a diagrammatic C/K-dual representation for a five-point 
building-block involving generalized fields. 

1.2. Color-kinematics duality for g*g*{s*s*) —>■ qqg 

As a non-trivial example of C/K-duality for dimension- 
ally regulated amplitudes, we consider again the process 
gg qqg, already discussed in Section 6 , but now re¬ 
garding the initial state gluons as d-dimensional particles 
(whereas the final state remains fully four-dimensional). 
This amplitude would contribute to a d-dimensional uni¬ 
tarity cut of a loop-level amplitude where the gluons pi 
and p 2 appear as virtual states. 

Within FDF, the full amplitude is obtained by combin¬ 
ing the contributions of three different processes involving 
generalized four-dimensional initial particles, 

9*9* -t 999: 



Figure 14: Jacobi combinations for FDF particles. 


s*s* -)> qqg, 

g's' -)> qqg. (94) 

However, thanks to the selection rules illustrated in Appendix A, 
g*s* —?► qqg vanishes and the problem decouples in the de¬ 
termination of the C/K-dual representation of two individ¬ 
ually gauge invariant amplitudes, 9 * 9 * qqg and s*s* —>■ 
qqg. The tree-level contributions to 9 * 9 * —>■ qqg are shown 
in Fig. 15. The Feynman diagrams for s*s* —>■ qqg can be 
easily obtained by replacing all generalized gluons lines g* 
with scalars s*. 

Since in both cases the number of graphs, the relations 
among their color factors and, as consequence, the set of 
constraints to be imposed on the shifted numerators (32) 
are exactly the same as in the example of Section 6 , we 
will simply discuss the relevant modifications to be taken 
into account in order to adapt the calculation to general¬ 
ized fields. 

We notice that, while for g*g* the redistribution of nu¬ 
merator ni 6 among cubic diagrams is still given by Eq. (47), 
for the s*s* case we have, as discussed in Section 7.1(b), 

ci6''t-16 = c^n^-iQ -b csngpe- (95) 


Nevertheless, for both processes the diagrammatic expan¬ 
sion of the amplitude can be still read from the r.h.s. of 
Eq. (49), provided the replacement 

sn -t (sn - M^) for i ^ 2 , 

S 2 i -t (s 2 * - P^) for i 1, (96) 

which account for internal massive propagators. These 
modifications affect the kinematic terms of the decomposi¬ 
tion (52), as well as the entries of the matrix (54). Finally, 
the anomalous terms that enter the definition of the vector 
(55) are given, for g*g* -)■ qqg, by 
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0502 \rA ' f 1 

OL^ 

gen 


•U4 eiders'" = S23I^[2,4.7], 


VAS'^^e'^^ — Si4V^[6,7,14]i 


^<>[1.3, 12 ] = ^3 [^q(Pl +P2,PA,Pd„Pb)a.a^ VAlV^^yVpaia^ (“Pi - P2,Pl,P2) sT + «3;16 = Sl2t/J[1.3,12], 

^<>[1.4.15] =U3 4'g'qq(P2,P4+P5,P3,Pl)„,„, (|>4 + l^s) ^5 
</'[2.4,7] = M 3 ^2 (I'a +1*2 + '^g'gq’q (P1:-P4,P3 +P 2 , 

</>[2.5.11] = '“3 [>^g*g*qq (pi +P5,P4,P3,P2)„„J ^4 (pi +P5,/^^) V/Jaias (“Pi - Pb,Pl,Pb) + M5;16 = Sl5</5[2,5,11] > 

<^[6,7.14] =M3 [-/g'gq'q (P5,P4+Pl,P3,P2)„,„, f]^4 + 1*1 + W 
</'[6.8,9] 

</^[9,10,15] 

</>[ll,10.13] =M 3 Jg*gq*q (Pl:-P4 +P2,P3,P5)a,a, (^ 4 +^ 9 +*^')'' 

^^]12,13,14] = 

and, for s*s* —>■ qqg, 


U 3 [4‘g‘qq (P2 +P5,P4,P3,Pl)aaJ M4n“f„ (P2 +P5,Ai^) V/Jasas [-P 2 - Pb,P 2 ,Pb) + M8;16 = S25</5]6,8.9], 

5^ 


M 3 

M 3 

M 3 


^1 Vr3 +1*1 '^g*gq*q(P5,P4,P3+Pl,P2)„,„, 




4*g*qq b2,P4,P3 +P5,Pl)c 


n4e2^e5® — si3V^]9,io,i5]i 
M4£i'^£^ 5'' = S24¥’]ll,10,13]! 


M4ei^£^2^ — S35(/9[12,13,14]) 


(97) 


<?i’]1.3.12] = M3 [Jq {pi +P2-,PA-,P3-,P5)aa^] VA^^ey (Pl - P-2)p = Sl2P[l,3,12], 


<P[lA,lb] = M 3 Js*s*qq (j>2,PA+P5,P3.,Pl) (^ 4 +^ 5 ) ^ 


1^4 = S45</5]1,4.15]; 


<(>]2.4.7] = M 3 7^ (^ 3+^2 + 


{Pl,PA,P3+P2,P3)c 


VaE^ — S23V5]2.4.7]i 


</^l2.5.11] — M3 [Js«s*qq (Pl + Pb,PA,P3,P2)] Va (2pi +P 5 ) ‘ ^5 + M5;16 — Sl5</^]2,5.11]; 


']6.7.14] — M 3 


<7s*gq*q 


{Pb,PA+Pl,P3,P2) 


«5 vr4+l^i+*^T'' 


-)A 


QIC 

M4e5 = S14</2]6,7,14] ) 


<(>]6.8.9] — M3 [Js«s*qq (P2 + PS, P4, P3, Pi)] M 4 (2^2 +P 5 ) ‘ ^5 + M8;16 — S25<P]6,8.9]; 


</']9,10,15] — M 3 
<(>]11.10.13] = M 3 


+ + *P7 


s*gq*q 


^s*gq*q (Pl,P4 +P2,P3,P5)c 


(P5,P4,P3 +P1,P2)„5 
+1>2 + *P7®) 7® 


M4 £ 5 ® — Sl3'P]9,10,15]i 
M4£5'* = S24<P]11,10,13]! 


<(>]12.13.14] =M 3 ^5 (^3 +^ 5 ) Js-s*qq(P2,P4,P3 +P5,Pl) 


M4 — S35<P]12.13.14]- 


(98) 


The systems AA = (p associated to both amplitudes still 
satisfy the condition (57), which ensures the existence of a 
C/K-dual representation, depending on four arbitrary pa¬ 
rameters, whose expression follows the structure of (68). 
As for the pure four-dimensional case, the analytic expres¬ 
sions of the dual numerators have been obtained for generic 
polarizations and numerical checks of the result have been 
performed for different helicity configurations, including 
longitudinal polarizations of generalized gluons. 

8. Conclusions 

In this letter we investigated, from a diagrammatic 
point of view, the off-shell color-kinematics duality for am¬ 
plitudes in gauge theories coupled with matter in four as 
well as d dimensions, within the Four Dimensional For¬ 
mulation variant of the Four Dimensional Helicity scheme. 


This duality, first observed at tree-level for on-shell four- 
point amplitudes, is non-trivially satisfied within higher- 
multiplicity tree-level or multi-loop graphs, due to pres¬ 
ence of contact terms which violate the Jacobi identity 
for numerators. We studied the source of such anomalous 
terms in gg -A ss,qq,gg scattering processes. Working in 
axial gauge, we have explicitly shown that, whenever the 
Lie structure constants obey a Jacobi identity, the anal¬ 
ogous combination of their kinematic numerators can al¬ 
ways be reduced to a sum of numerators of sub-diagrams, 
with one or two denominators less. 

Our decomposition provides a systematic classification 
of the duality-violating terms into a reduced number of ef¬ 
fective vertices and, since they vanish when on-shellness is 
imposed on the four particles identifying the Jacobi rela¬ 
tion, it immediately allows to recover the color-kinematics 
duality for multi-loop cut-integrands built from Feynman 
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Figure 15: Feynman diagrams for g*g* —)• qqg. The contributions to 
s* s* qqg are obtained by replacing all generalized gluons g* with 
s* lines. 

rules. 

We consider this study as an independent step towards 
a different, yet direct perspective to the diagrammatic un¬ 
derstanding of color-dual graphs. Our approach, based 
on the direct inspection of Feynman diagrams and the 
identification of a set of constraints able to remove C/K- 
violating terms, offers a method for the construction of 
dual numerators which is alternative to the traditional one, 
where, starting from general ansatz on the functional de¬ 
pendence of numerators on external momenta and polar¬ 
izations, Jacobi-like symmetries are imposed. 

Off-shell recurrence relations, in tandem with general¬ 
ized gauge transformations, can play an important role 
for gaining further understanding of the on-shell C/K- 
duality in higher-multiplicity processes, as shown by our 
explicit determination of dual numerators for the tree-level 
gg —>■ qqg amplitude, first computed in four dimensions 
and later in d dimensions, where the initial state gluons 
were considered as dimensionally regulated particles. 

We expect this approach to have a natural extension 
at loop-level, which will be object of future studies. 
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Appendix A. 

In this Appendix we briefly recall the main features of 
the FDF scheme and we give the set of Feynman rules used 
in the calculations discussed in Section 7. 

• We use barred notation for quantities referred to un¬ 
observed particles, living in a d-dimensional space. 
Thus, the metric tensor 

r’'= 9^'^+r'', (A.i) 

can be decomposed in terms of a four-dimensional 
tensor g and a — 2e-dimensional one, g. The ten¬ 
sors g and g project a d-dimensional vector q into 
the four-dimensional and the — 2 e-dimensional sub¬ 
spaces, respectively. 

• d-dimensional momenta i are decomposed as 

i = e + e, F = = (a. 2 ) 

• The algebra of matrices 7 ^ = g^^ 7 '", 

[ 7 “, 7 ^] = 0, { 7 “, 7 ''} = 0, (A.3a) 

= (A.3b) 

is implemented through the substitutions 

(A.4) 

together with the set of selection rules, (—2e)-SRs, 


qABqBC 

= G^^, 

qAA 

= 0 , 

qAB _ 

pA ^AB 

= F^, 

^A^^A 

= 0 , 

= 1 , 




= 1 . 

(A.5) 


which, ensuring the exclusion of the terms containing 
odd powers of g, completely defines the FDF and 
allows the construction of integrands which, upon 
integration, yield to the same result as in the FDH 
scheme. 

• The spinors of a d-dimensional fermion fulfil the com¬ 
pleteness relations 

u\ (i) u\ (£) = /-!- f/ry® -I- m , 

A=± 
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(A.6) 


V\ {i) vx (£) = / + - m , 

\=± 


which consistently reconstruct the numerator of the 
cut propagator. 

• In the axial gauge, the helicity sum of a d-dimensional 
transverse polarization vector can be disentangled in 


d-2 



i=l 


{£, v) £i{d) (£, n) = 




(A. 7 ) 


where the first term can be regarded as the cut 
propagator of a massive vector boson, 


E = + — , (A.8) 

A=±.0 ^ 


whose polarizations obey the expected properties 


4 (£)= 0 , 
£2(£) = -1, 
ex(£)-£= 0 


e±(£) ■ ezp(£) — —1, 

e±(£)-so(£)= 0, 

(A.9) 


For the explicit expression of polarization vectors as 
well as generalized spinor we refer the reader to [54] . 
The second term of the r.h.s. of Eq. (A.7) is related 
to the numerator of cut propagator of the scalar s* 
and can be expressed in terms of the (—2e)-SRs as: 

+ G^^-Q^Q^. (A.IO) 


Within FDF scheme, the QCD d-dimensional Feynman 
rules in axial gauge have the following four-dimensional 
formulation: 


k 

•JLSLOJLOJU 

b,!3 


gab 

— IJ? + iO 





(A.lla) 


gab 


,A b, B 


— g? + iO 


(pAB _ qAqB^ ^ 


= i 


^ + i/iy® -|- TO 

k2 _ jyg2 _ ^2 _|_ jQ 


(A.llb) 

(A.llc) 



= - 3 /“'“ [{ki-k2Vg<^P 

+ (fc2-fc3)V^ + (A:3-A:i)V“] , 

(A.lld) 


2, b, B 

=- 3 r'’ 4 fc 2 -fc 3 )“C?^‘^, (A.lle) 

3, c, C '' 




(fci =0, h = ±£). (A.Ilf) 

= -ig^ r'"’" 

_|_ ja;ac jxbd _ g°‘^ 

(A.llg) 

_ ^^2 gOiS ^jxab jxcd _|_ j^xac j^xbd'^ 



(A.llh) 

(A.lli) 



= -*d (i% 4r- 


(A.llj) 


In the Feynman rules (A.11) all the momenta are incom¬ 
ing and the scalar particle s* can circulate in the loop 
only. The term appearing in the propagators (A.lla)- 
(A.llc) enter only if the corresponding momentum k is 
d-dimensional, i.e. only if the corresponding particle cir¬ 
culates in the loop. In the vertex (A.Ilf) the momentum ki 
is four-dimensional, while the other two are d-dimensional. 
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